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Abstract 

In this work the possibihty of detecting solar electron antineutrinos produced by 
a solar core magnetic field from the KamLAND recent observations is investigated. 
We find a scaling of the antineutrino probability with respect to the magnetic field 
profile in the sense that the same probability function can be reproduced by any 
profile with a suitable peak field value. In this way the solar electron antineutrino 
spectrum can be unambiguosly predicted. We use this scaling and the negative results 
indicated by the KamLAND experiment to obtain upper bounds on the solar electron 
antineutrino flux. We get (pp < 3.8 x 10-^(j){^B) at 95% CL. For 90% CL this becomes 
(pp < 3.4 X 10~^(/>(*i?), an improvement by a factor of 3-5 with respect to existing 
bounds. These limits are independent of the detailed structure of the magnetic field 
in the solar interior. We also derive upper bounds on the peak field value which are 
uniquely determined for a fixed solar field profile. In the most efficient antineutrino 
producing case, we get (95% CL) an upper limit on the product of the neutrino 
magnetic moment by the solar field fiB < 2.8 x 10~^^ MeV or Bq < 4.9 x lO^G for 
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1 Introduction 



The recent results from the KamLAND experiment [1] have asserted that the large mixing 
angle solution (LMA) is the dominant one for the 34 year old solar neutrino problem SNP 
[2]. Although neutrinos were known, before KamLAND data, to oscillate [3, 4], it was not 
clear if neutrino oscillations were the major effect underlying the solar neutrino deficit or 
whether they played any role at all. It had been clear however that this deficit had to rely 
on 'non-standard' neutrino properties. To this end, the spin flavor precession (SFP) [5, 6, 7], 
based on the interaction of the neutrino magnetic moment with the solar magnetic field was, 
second to oscillations, the most attractive scenario [8]. 

SFP, although certainly not playing the major role in the solar neutrino deficit, may still 
be present as a subdominant process, provided neutrinos have a sizeable transition magnetic 
moment. Its signature will be the appearance of solar antineutrinos [6, 9, 10] which result 
from the combined effect of the vacuum mixing angle 9 and the transition magnetic moment 
liy converting neutrinos into antineutrinos of a different fiavor. This can be schematically 
shown as 

VeL J^HL ^en, (1) 
T^eL '^MB ^ (2) 

with oscillations acting first and SFP second in sequence (1) and in reverse order in sequence 
(2). Oscillations and SFP can either take place in the same spatial region, or be spatially 
separated. Independently of their origin, antineutrinos with energies above 1.8 MeV can be 
detected in KamLAND via the observation of positrons from the inverse /3-decay reaction 
+ P ^ + e"^ and must all be originated from neutrinos. 

The purpose of this work is to relate the solar magnetic field profile to the solar an- 
tineutrino event rate in KamLAND which is a component of the total positron event rate in 
the reaction above. In a previous paper [11] the question of what can be learned about the 
strength and coordinate dependence of the solar magnetic field in relation to the current 
upper limits on the solar z/g fiux was addressed. The system of equations describing neu- 
trino evolution in the sun was solved analytically in perturbation theory for small fi^B, the 
product of the neutrino magnetic moment by the solar field. The three oscillation scenarios 
with the best fits were considered, namely LMA, LOW and vacuum solutions. In particular 
for LMA it was found that the antineutrino probability depends only on the magnitude of 
the magnetic field in the neutrino production zone. Neutrinos were, in the approximation 
used, considered to be all produced at the same point {x — 0.05Rs), where neutrino 
production is peaked. In this work we will consider the more realistic case of a convolution 
of the production distribution spectrum with the field profile in that region. It will be 
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seen that this convolution leads to an insensitiveness of the antineutrino probability with 
respect to the solar magnetic field profile, in the sense that different profiles can correspond 
to the same probability function, provided the peak field values are conveniently scaled. 
As a consequence, an upper bound on the solar antineutrino fiux can be derived which is 
independent of the field profile and the energy spectrum of this fiux will also be seen to be 
profile independent. 

Up to now, in their published data from a 145 day run, the KamLAND [1] experiment 
has observed a total antineutrino flux compatible with the expectations coming from the 
nearby nuclear reactors. Given this fact, and evaluating the positron event rate for the 
above reaction for different solar fleld proflles, we will derive upper bounds for the peak 
fleld value in each profile. In our study we will assume the astrophysical upper bounds on 
the neutrino magnetic moment < (1 — 3) x 10~^^//b [12] to be all satisfied and take 
= 10-^Vs- 

Upper limits on the solar antineutrino fiux, the intrinsic magnetic moment and the 
magnetic field at the bottom of the convective zone were recently obtained [10] from the 
published KamLAND data. Here we address however a different antineutrino production 
model where the magnetic fleld at the solar core is the relevant one. 



2 The solar antineutrino probability 

We start with the probability that a Ugj^ produced inside the sun will reach the earth as a 



P{Ue^ Vej,) = P{Ue^ ^ Z/^^; Rs) X P{u^^ ^ Z/^^; Res) (3) 



in which the first term is the SFP probability, Rg is the solar radius and the second term is 
given by the well known formula for vacuum oscillations 

P{^^^u ^ ^e^; Res) = sm'2e sm' {-j^Resj = 2 ■ (4) 

Here R^s is the distance between the sun and the earth and the rest of the notation is 
standard. Since l.SMeV <E < IbMeV and, for LMA, Am^ = 6.9 x lO'^eV^, sin^ 29 = 1 
[1] , we take the P^^ — > Pea vacuum oscillations to be in the averaging regime. 
The SFP amplitude in perturbation theory for small /iB is [11] ^ 

^'"^ ^ = — — ■ 



^For notation we refer the reader to ref. [11]. 
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A key observation is that the antineutrino appearance probabihty is dependent on the 
production point of its parent neutrino so that the overall antineutrino probability is 

P{Ue^ Uen) = \ j IM^ei^ ^ I^hr)? f B{ri)dri (6) 

where Jb represents the neutrino production distribution function for Boron neutrinos [13] 
and the integral extends over the whole production region. As shall be seen, owing to this 
integration, the energy shape of probability (6) is largely insensitive to the magnetic field 
profile. 

The positron event rate in the KamLAND experiment originated from solar antineutri- 
nos is then 

/•oo i-Em , , 

S = Qo dEe e{EjR{Ee, EJME)a{E)dE. (7) 

JEO jEm 

In this expression Qq is a normalization factor which takes into account the number of atoms 
of the detector and its live time exposure [1] and E is the antineutrino energy, related to 
the physical positron energy hy E'^ = E — {mi^ — mp) to zero order in 1/M, the nucleon 
mass. We thus have Em = l.SOAMeV, while the KamLAND energy cut is E° = 2.6MeV. 
The functions e and R denote the detector efficiency and the Gaussian energy resolution 
function of the detector 



RiE„E^)^^=exp 

SVZTT 



2s2 



(8) 



In our analysis we use for the energy resolution in the prompt positron detection the expres- 
sion s{Ee) = 0.0062 + 0.065y'i?e with all energy units in MeV. This is obtained from the raw 
calibration data presented in Ref.[14]. Moreover, we assume a 408 ton fiducial mass and the 
detection efficiency is taken independent of the energy [14], e ~ 80%, which amounts to 162 
ton. yr of antineutrino data. The antineutrino cross section a{E) was taken from ref.[15] 
and we considered energy bins of size E^, = 0.425 MeV in the KamLAND observation range 
(2.6 — 8.125) MeV [1]. The antineutrino spectral flux 4>u{E) in eq. 7 can be written as 
(j)p{E) — (fp- y. f{E) where 0° is the total antineutrino flux and f{E) is some function of 
the energy normalized to one. The function 4>v{E) is on the other hand a simple product 
of the Boron neutrino spectral flux (psiE) which can be found in Bahcall's homepage [13] 
and the antineutrino appearance probability we obtained above: (f>p{E) = (pB^E) x P{E). 
The almost insensitivity of the shape of P{E) to the shape of the magnetic field profile is 
thus necessarily refiected in (f)p{E). The only significant dependence appears on the normal- 
ization constant 0° which is essentially proportional to the square of the magnetic field at 
the solar core. We make use of this behavior to obtain, for each given profile, upper limits 
on the core magnetic field, the total antineutrino flux and the intrinsic neutrino magnetic 
moment. 
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As mentioned above, for the LMA solution only the solar field profile in the neutrino 
production region [11] can affect the antineutrino flux. Hence we will discuss three proflles 
which span a whole spectrum of possibilities at this region. We study from a vanishing fleld 
(proflle 1) to a maximum fleld at the solar center, with, in this second case, either a fast 
decreasing fleld intensity (proflle 2) or a nearly flat one (proflle 3) in the solar core (see flg. 
1, lower panel). Thus, we consider respectively the following three profiles 

Profile 1 

B{r) = So [cosh (9r) - 1] , |r| < Tc (9) 

B{r) = 5o/cosh[25(r - tr)] , |r| > r^, (10) 

with Tc = 0.08, tr = 0.16, 
Profile 2 

B(r) = So/cosh(15r) , |r| > 0, (11) 

Profile 3 

B{r)^Bo[l-{r/r,f] , [r| < (12) 

with Tc = 0.713. 

We also show in fig. 1 (upper panel) the production distribution spectrum, so that 
a comparison between the strength of the field and the production intensity can be directly 
made. 

The antineutrino production probabilities as a function of energy for each of these 
proflles are given in flg. 2. In the flrst panel, the values of the peak fleld are chosen so 
as to produce a flxed number of events. In this case the probability curves differ only 
slightly in their shapes while their normalizations are the same. The curves are in any 
case similar to the SFP survival probability ones [16] in the same energy range. In the 
second panel of fig. 2 the antineutrino probabilities for a common value of the peak field 
and these three different profiles are shown. It is hence apparent from these two graphs 
how the distribution of the magnetic field intensity is determinant for the magnitude of the 
antineutrino probability, but not for its shape. One important reason for this behavior is 
that we have integrated the antineutrino probability over the Boron production region. 



3 Results and discussion 

The antineutrino signal for any magnetic field profile B{r) can be written, taking into 
account the previous formulas and the near invariance of the probability shape (see fig. 2), 
as 

S^[B{r)] = aS'^ (13) 
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where ^S"^ is the antineutrino signal taken at some nominal reference value Bq for the field at 
the solar core for a certain reference profile B^. This profile dependent parameter a, being 
a ratio of two event rates given by eq.(7) for different profiles, can thus be simplified to 




where the integrals extend over the production region. As we mentioned before, for con- 
creteness we have fixed along this discussion the neutrino magnetic moment /ii, — 10~^^fj,B- 

We will now obtain bounds on parameter a and the peak field Bq for each profile derived 
from KamLAND data, applying Gaussian probabilistic considerations to the global rate in 
the whole energy range, Ej, — (2.6 — 8.125) MeV, and Poissonian considerations to the event 
content in the highest energy bins {Ee > 6 MeV) where KamLAND observes zero events. 
We denote by Sp the event rate with Bq = lO'^G for each given profile {Sp = S'p(lO^G)). 
Taking the number of observed events and subtracting the number of events expected from 
the best-fit oscillation solution [{Am'^ , sin^ 29) lm a = (6.9 x 10^^ ey^, 1)] and interpreting 
this difference as a hypothetical signal coming from solar antineutrinos, we have 

Sr = S,,,-Sreact{LMA). (15) 

Inserting [14] Sobs = 54.3±7.5 and Sreact{LMA) = 49± 1.3, we obtain Sobs - Sreact = aS^ < 
17.8 (20.2) at 90 (95)% CL. Within each specific profile it is seen from (14) that the quantity 
a is simplified to o; = (Sq/IO^G)^, so that the previous inequality becomes 

Bl < ^(lO^G)'- (16) 

In this way we can derive for each given profile an upper bound on Bq. The quantity S^ for 
profiles 1, 2 and 3 and the respective upper bounds on Bq are shown in table 1. These upper 
limits can be cast in a more general way if do not fix the neutrino magnetic moment. To this 
end we will consider an arbitrary reference value = IQ^^'^^b- Then within each profile, 
a = {fi,yBo/fi^ lO^G)^, where in the numerator and denominator we have respectively the 
peak field value and some reference peak field value of the same profile. In the same manner 
as before we can derive the upper bounds on h^Bq which are also shown in table 1. 



From the definition of a (14) it follows that the upper bounds on the antineutrino 
flux are independent of the field profile. These turn out to be 0p < 0.00340(^5) and 
< 0.00380(^5) for 90 and 95% CL respectively. 

We can similarly and independently apply Poisson statistics to the five highest energy 
bins of the KamLAND experiment. No events are observed in this region and the expected 
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signal from oscillating neutrinos with LMA parameters is negligibly small. We use the fact 
that the sum of Poisson variables of mean fii is itself a Poisson variable of mean J2 A^i- The 
background (here the reactor antineutrinos) and the signal (the solar antineutrinos) are 
assumed to be independent Poisson random variables with known means. If no events are 
observed and in particular no background is observed, the unified intervals [17, 18] [0, ecL] 
are [0, 2.44] at 90% CL and [0, 3.09] at 95% CL. 

Prom here, we obtain aS^ < ecL or a < ecL/Sp. Hence, as in the previous case, we 
have 

B'o < 'j^ilO'Cf. (17) 

Using the expected number of events in the first 145 days of data taking and in this energy 
range (6 — 8.125) MeV, we have derived upper bounds on Bq (90 and 95% CL) for all 
three profiles. They are shown in table 2 along with the upper bounds on /IuBq taking 
III, as a free parameter. The antineutrino flux upper bounds are now 0p < 0.00490(^5) 
0p < 0.00550(^5) at 90 and 95% CL respectively. The KamLAND expected signal for an 
arbitrary held proflle corresponding to 95% CL is shown in flg. 3. 

The differences in magnitude among the bounds on Bq and ij,,^Bq presented in tables 1 
and 2 for the different profiles are easy to understand. In fact, recalling that the ^B produc- 
tion zone peaks at 5% of the solar radius and becomes negligible at approximately 15% (fig. 
1), then in order to generate a sizeable antineutrino fiux, the magnetic field intensity should 
lie relatively close to its maximum in the range where the neutrino production is peaked. 
Thus for profile 1 the value of Bq required to produce the same signal is considerably larger 
than for the other two, while profile 3 is the most efficient one for antineutrino production. 

As referred to above, for different field profiles the probability curves will differ only 
slightly in their shape if they lead to the same number of events. In other words, for a given 
number of events the probability curves are essentially the same, regardless of the field 
profile, a fact illustrated in fig. 2. As a consequence, the energy spectrum of the expected 
solar antineutrino flux will be nearly the same for any profile. In fig. 4 we plot this profile 
independent spectrum together with the ^B one [13], so that a comparison can be made 
showing the shift in the peak and the distortion introduced. 

4 Conclusions 

To conclude, now that SEP is ruled out as a dominant effect for the solar neutrino deficit, 
it is important to investigate its still remaining possible signature in the solar neutrino 
signal, namely an observable z/g flux. Our main conclusion is that, from the antineutrino 
production model expound here, an upper bound on the solar antineutrino flux can be 
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derived, namely < 3.8 x 10^^(t){^B) and 0,? < 5.5 x lO"^0(^i?) at 95% CL, assuming 
respectively Gaussian or Poissonian statistics. For 90% CL we found 0p < 3.4 x 10^^0(^-8) 
and 0p < 4.9 x 10^^0(^-8) which shows an improvement relative to previously existing 
bounds from LSD [19] by a factor of 3-5. These are independent of the detailed magnetic 
field profile in the core and radiative zone and the energy spectrum of this flux is also 
found to be profile independent. We also derive upper bounds on the peak field value which 
are uniquely determined for a fixed solar field profile. In the most efficient antineutrino 
producing case (profile 3), we get (95% CL) an upper limit on the product of the neutrino 
magnetic moment by the solar field n^B < 2.8 x 10"^^ MeV or Bq < 4.9 x lO'^'G for 
— IO^^^Hb- a recent study of the magnetic field in the radiative zone of the sun has 
provided upper bounds of (3-7) MG [20] in that region in the vicinity of 0.2 Rs which are 
independent of any neutrino magnetic moment. Therefore we can use them in conjunction 
with our results to obtain a limit on /x,^. Using i?o ~ 3 — 7MG, we get from the results for 
profiles 1-3: h < 0.7 — 9.6 x 10~^^/xs. Moreover, from the limits obtained in this work, if 
the 'true' solar profile resembles either a profile like 1 or 3, this criterion imphes that SFP 
cannot be experimentally traced in the next few years, since the peak field value must be 
substantially reduced in order to comply with this upper bound, thus leading to a much too 
small antineutrino probability to provide an observable event rate ^. On the other hand, for 
a profile like 2 or in general any one resembling a dipole field, SFP could possibly be visible. 
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Profile 




5o(90%CL) 


5o(95%CL) 


fi,Bo{90%CL) 


fi,Bo{95%CL) 






G 


a 


MeV 


MeV 


1. 


0.006 


5.27 X 10^ 


5.62 X 10^ 


3.05 X 10'^^ 


3.25 X 10-^^ 


2. 


0.137 


1.14 X 10^ 


1.21 X 10*^ 


6.60 X 10^19 


7.04 X 10-19 


3. 


0.224 


8.92 X 10^ 


9.50 X 10^ 


5.16 X 10-19 


5.50 X 10-19 



Table 1: Solar antineutrino event rates, upper bounds on the peak field value for jjiy = 
10-i^/iB and on ix^Bq for arbitrary and Bq, assuming Gaussian statistics in the whole 
KamLAND spectrum. 



Profile 


Sl{WG) 


5o(90%CL) 


5o(95%CL) 


/x,5o(90%CL) 


/x,5o(95%CL) 






G 


G 


MeV 


MeV 


1. 


0.004 


2.53 X 10« 


2.85 X 10^ 


1.47 X 10-1^ 


1.65 X 10-1^ 


2. 


0.079 


5.56 X 10^ 


6.25 X 10^ 


3.22 X 10-19 


3.62 X 10-19 


3. 


0.130 


4.34 X 10^ 


4.88 X 10^ 


2.51 X 10-19 


2.82 X 10-19 



Table 2: Same as table 1 assuming Poissonian statitics in the KamLAND energy range 
Ee = (6-8.125) MeV. 
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Figure 1: Upper panel: neutrino production spectrum (in arbitrary units) as a function 
of the radial coordinate. Lower panel: the three solar field profiles considered in the main 
text normalized to Bq, the peak field value. 
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Figure 2: Antineutrino probabilities for solar field profiles 1, 2 and 3. Upper panel: the peak 
field is chosen in each case so as to produce the same event rate in KamLAND, (see the 
main text). Lower panel: the same value of the peak field (Bq = lO^G) is seen in each case 
to lead to probabilities of quite different magnitudes. 
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Figure 3: The solid squares represent the MC expectation of the KamLAND positron spec- 
trum from reactor antineutrinos with no oscillations and the points with error bars represent 
the measured spectrum (from Fig. 5 in Ref.[l]). Solid triangles represent the positron spec- 
trum from solar antineutrinos (multiplied by 5) assuming profile 3 with peak field given by 
its 95% CL upper limit (Bq = 4.88 x lO'^G). All curves refer to the same time exposure of 
145 days. 
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Figure 4: The expected solar antineutrino spectrum and the neutrino one [13], both 
normalized to unity, showing the peak shift and the distortion introduced by the antineutrino 
probability. 
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